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Abstract 

^ ■ ^a(I) symmetry breaking is studied by introducing the flavor mixing interaction pro- 

posed by Kobayashi, Maskawa and 't Hooft. Combining the one gluon exchange inter- 
action, the rainbow like Schwinger-Dyson equation and the ladder like Bethe-Salpeter 
equation are derived. The anomalous PC AC relation in the framework of this approxima- 
tion is considered. The masses of the pseudoscalar mesons tt, r\ and rj' are calculated. It is 
found that the pion mass is not sensitive to the strength of the flavor mixing interaction. 
On the other hand, the masses of n and n' are reproduced by a relatively weak flavor 
mixing interaction, for which the chiral symmetry breaking is dominantly induced by the 
soft-gluon exchange interaction. The decay constants are calculated and the anomalous 
PCAC relation is numerically checked. It is found that the flavor structures of the n and 
n' mesons significantly depend on their masses and therefore it is quetionable to define a 
flavor mixing angle for rj and n' ■ 
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1 Introduction 



It is known that the classical QCD lagrangian is invariant under the C/l(3) x Ur(3) symmetry 
except for the quark mass term, and this symmetry is broken down to the LV(3) spontaneously 
in the low-energy QCD. In this case, the number of the Nambu-Goldstone bosons should be 
nine. However, the number of the observed light pseudoscalar mesons is eight. The ninth 
pseudoscalar meson, rf meson, is heavier than the other octet pseudoscalar mesons (tt°, tt + , tt~ , 
K + ,K-,K°, K and rj) which are well identified with the Nambu-Goldstone bosons. This is 
the well-known {7a (1) problem It is solved by realizing that the {7a (1) symmetry is broken 
by the anomaly. The phenomena related to the £7a(1) anomaly in the low-energy QCD have 
been studied in the following approaches. 

The first one is the 1/Nc expansion approach 0, ||. The key point is that the effect of the 
{7a(1) anomaly is higher order in the 1/Nc expansion and the low-energy effective lagrangian 
of QCD has been derived Recently, the expansion in powers of 1/Nc, momenta and quark 
masses was extended to the first non-leading order || and the reasonable description of the 
nonet pseudoscalar mesons was obtained. 

The second is the instanton approach ||. The instanton is a classical solution of the Eu- 
clidean Yang-Mills equation and may contribute a large weight in the Feynman path integra- 
tion. In the presence of the light quarks, instantons are associated with fermionic zero modes 
which give rise to the {7a (1) symmetry breaking. In the dilute instanton gas approximation, 
the {7a(1) breaking 6-quark flavor determinant interaction was derived in the three flavor case 
0. This approach has been developed to the instanton liquid picture of the QCD vacuum |J. 
In this picture, the instanton plays a crucial role in understanding not only the {7a(1) anomaly 
but also the spontaneous breaking of the chiral symmetry itself. 

In the third approach, the effective low-energy quark models of QCD were used to study the 
structure of the hadrons. The introduction of the instanton induced six-quark interaction to 
the effective quark model is one of the handy way to incorporate the the {7a (1) breaking effects 
in the low-energy effective quark model of QCD. The Nambu-Jona-Lasinio (NJL) model |§ is 
one of the simplest and widely used model in studying the structure of the Nambu-Goldstone 
bosons. Using the three-flavor NJL model with the instanton induced six quark interaction, 



2 



properties of the nonet pseudoscalar mesons were investigated [Hj. A shortcoming of this 
approach is that the rj' mass has unphysical imaginary part associated with the unphysical 
decay channel rf — > qq. Recent study of the radiative decays of the r\ meson in the NJL model 
has shown that the observed values of the mass and radiative decay amplitudes are reproduced 
well with a rather strong Ua(1) breaking interaction pTTJ . Such strong {7a(1) breaking may be 
consistent with the instanton liquid picture of the QCD vacuum. 

In contrast with the instanton liquid model, the study of the QCD Schwinger-Dyson (SD) 
equation for the quark propagator in the improved ladder approximation (ILA) has shown 
that the spontaneous breaking of the chiral symmetry is explained by simply extrapolating the 
running coupling constant from the perturbative high-energy region to the low-energy region 

Then, the Bethe-Salpeter (BS) equation for the J PC = 0~ + qq channel has been solved in 



the ILA and the existence of the Nambu-Goldstone pion has been confirmed ||13|| . The numerical 
predictions of the pion decay constant and the quark condensate (ipip) are rather good. It 
has been also shown that the BS amplitude shows the correct asymptotic behavior as predicted 
by the operator product expansion (OPE) in QCD ]14|. The masses and decay constants for 
the lowest lying scalar, vector and axial-vector mesons have been evaluated by calculating the 
two point correlation functions for the composite operators ipMip. The obtained values are in 



good agreement with the observed ones [15]. 

Recently, the current quark mass term has been introduced in the studies of the BS ampli- 



tudes in the ILA JT6fl and the reasonable values of the pion mass, the pion decay constant and 
the quark condensate have been obtained with a rather large Aqcd- It has been also shown 
that the pion mass square and the pion decay constant are almost proportional to the current 
quark mass up to the strange quark mass region. 

Since the rj and rj' system is expected to be sensitive to the Ua{X) anomaly, the study of the 
r\ and rj' structure may give us information on the roles of the anomaly in the low-energy 

QCD. The purpose of this paper is to study the properties of the r] and rj' mesons by solving 
the coupled channel BS equation in the ILA. The effect of the {7a(1) anomaly is introduced by 
the instanton induced six-quark determinant interaction. The instanton size effects are taken 
into account by the form factor of the interaction vertices. It guarantees the right asymptotic 
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behavior of the solutions of the SD and BS equations. 

There have been many studies of the pion BS amplitude using the effective models of QCD 
and /or the approximation schemes of QCD ||17|| . As for the r\ and rj' system, Jain and Munczek 
model |18|] has been applied to them ||19fl . They have introduced the effect of the £7a(1) anomaly 



by simply adding the additional mass term in the flavor singlet pseudoscalar meson channel by 
hand and the reasonable values of the masses and decay constants have been obtained. 

It is known that the introduction of the two gluon exchange diagrams in the calculation of the 
f] and 7]' BS amplitudes beyond the ladder approximation do not break the Ua(1) symmetry with 
the perturbative gluon propagator. Recently, it has been shown that if the gluon propagator 
has the strong infrared singularity, the £/a(1) symmetry breaks PC|| . The relaton between this 
approach and the instanton approach is not clear. 

The paper is organized as follows. In Sec. |2| we explain the model Lagrangian we have 
used in the present study. In Sec. |]the Conwall-Jackiw-Tomboulis (CJT) effective action pl | 
calculated from our model Lagrangian is presented. The SD and BS equations are derived from 
the CJT effective action in Sec. f|and Sec. |5|. In Sec. ||the meson decay constant is derived. In 
Sec. |7| the Nambu-Goldstone solution of the BS equation is derived and the anomalous PCAC 
relation is discussed in Sec. ^|. Sec. [5] is devoted to the numerical results. Finally, summary and 



concluding remarks are given in Sec. 10 



2 Model 

We use the flavor three (Np = 3) effective model whose lagrangian density is given by 

:=^f(d 2 )(i@-m )^ + £ GE ty^]+£ F M&^}, (1) 
i;:=(u,d 7 s) T (2) 
where /(£) is a cut-off function defined by [l^ 



/(£) = 1 + M0(£-A2 V ), M^oo. (3) 

m denotes a diagonal quark mass matrix m = diag(m g , m g , m s ) under the assumption of the 
isospin invariance. £qe denotes a gluon exchange interaction 

I Jpp'qq 1 
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where / p denotes / and represents the Euclidean momentum. The indeces m, n, ... 
represent combined indeces in the color, flavor and Dirac spaces. In Eq.([|) we employ the 
Landau gauge gluon propagator, 

k»k v \ -1 



iD"(k) = ^ - //2 ; //2 . 
and the Higashijima-Miransky type running coupling constant g 2 defined as follows. 



(6) 



9\P% <&) = 0(p 2 E ~ q E )g\p E ) + 9(q E ~ P E )g\q E 



(7) 



with 



9 2 (pl 



1 1 



2(3 (1 + tip) 5 



3t 



IF — CO 



tlF — to 



1 3tip — to -\- 2, 
{ Wo (1+tiF) 2 



for t IF < t 



for to <t < tip 



for t < t 



iv QCD 



1 11JV C - 2iV 



(9) 



(10) 



(4tt) 2 3 

In Eq.(||) the infrared cut-off tip is introduced. Above tip, g 2 {p E ) develops according to the one- 
loop result of the QCD renormalization group equation and below t , g 2 (p%) is kept constant. 
These two regions are connected by the quadratic polynomial so that g 2 {p E ) becomes a smooth 
function. Here Nq is the number of colors and Nf is the number of active flavors. We use 
Nc = Nf = 3 in our numerical studies. 
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£fm is the £7a(1) symmetry breaking flavor mixing interaction (FMI), our interest, given by 
Cn&tyMx) = \G D e hhh e^w{d Xx -d yi -d X2 -d V2 -d X3 -d y3 ) 



(11) 



x { Wg^MhiVi)] ^(^2)^/2(2/2)] ^3(^3)^/3 (ya)] 

+ 3[V» S1 {Xl)ll>h (Vl)) bPg 2 (X 2 )l5^f 2 ^2)] [iPgz (^3)75^/3 (s/a)] 



where /1, #1, • • • are flavor indices, e denotes the antisymmetric tensor with e uds = 1 and the 
symbol * at the end of the equation means x±,yi, ■ ■ ■ — > x after all derivatives are operated. 
This type of the £7,4 (1) symmetry breaking six-quark interaction has been introduced in Ref. 



22] before the discovery of the instanton induced interacion m. There is a minor difference 
of the flavor-spin structure between £fm and 't Hooft instanton induced interaction. In the 
rainbow-ladder approximation, these two interactions give the same effects on the rj-r)' system. 

We introduce a weight function w(- ■ ■) which is necessary so that FMI is turned off at the 
high energy. We use the following separable Gaussian form 

/02 + Q2 + Q2 + Q2 + Q2 + Q2 N 

w(d Xl ;d yi ;d X2 ;d y2 ;d X3 ;d m )=w[ x ^ ^ + x > + , (12) 

w(fj, 2 ) := exp(— k/j 2 ). (13) 

This weight function is convenient for a numerical calculation as it satisfies the association rule 

w(-p 2 - q 2 - k 2 ) = w(-p 2 )w(-q 2 )w(-k 2 ). (14) 

But this particular set of the momenta in the argument of the weight function modifies the 
form of the Noether current for the axial-vector transformation. This is the same problem 
occurred in the Higashijima-Miransky approximation in the £ge term discussed extensively in 



Ref.|16[. The explicit form of the Noether axial- vector current in this model is very complicated 
and we do not show it. We treat the exact Noether axial- vector current within the ladder like 
approximation. We will show that the (modified) Ward-Takahashi identity for axial-vector 



current and the PCAC relation holds. This approach is studied in Ref.f23f 



On the other hand, if one does not want to modify the Noether current, one has to employ 
an appropriate form of the argument of the running coupling constant, such as 

? { ... )= ?l(?E±l!k-«E±3k) 2 \ (15 ) 




in £qe and similarly 

w{- ..) = w ^ dxi + 9x2 + 8x ' i 2 ~J ) Vl ~ ^ ~ 9m 1 I (16) 
in £fm- 

In our model, there are nine axial- vector currents, J^(a = 0, •••,8), which satisfy the 
anomalous PC AC relation 

d»J^(x) = 2[m J 5 ] a (x) + S a0 A(x), (17) 

[m J 5 \ := ^Z7 5 V (18) 

where X a denotes the Gell-Mann matrix in the flavor space. A(x) corresponds to the explicit 
C/a(1) symmetry breaking and is proportional to G F in Eq. (O) . In QCD, A(x) is given by 
3a s 



Sty 



e ,upa F a F a (lg) 



3 Effective Action 

To derive the Schwinger-Dyson equation and the Bethe-Salpeter equation, we use the Cornwall- 
Jackiw-Tomboulis (CJT) effective action formulation JH]]. In Ref.[IB|, we have already derived 
the CJT effective action in the lowest order (rainbow-ladder) approximation in the framework 
of the ILA model. Here we add a new term I^m^f] which contains the lowest order effect of 
the flavor mixing interaction (FMI). 

T[S F ] := iTrLn[S F ] - iTr^Sp] + T GE [S F ] + T FM [S F ] (20) 

Tce^f] corresponds to the two-loop (eyeglass) diagram using gluon exchange interaction and 



is defined by Eq.(15) in Ref. [16|. Since FMI is a six-quark interaction, the lowest two particle 
irreducible vacuum diagram of FMI is a three loop (clover) diagram. For simplicity, we take 
only the dominant term in 1/Nc expansion for Tpm^f] term as 

T fm [Sf] = ^ / d'xe^e^w ^% +%+ + ^ (21) 
x | - tr< DC ) [S Ffigi ( Vl , Xl )}trW { s Fhg2 (y 2} X2 )}trW [S Ffaga (y 3 ,x 3 )] 
- 3tr( DC )[S F/l9l ( yi ,xO]tr( D % 5 ^^^^ 




Figure 1: The diagram for the SD equation. 

In this approximation, the global SUl(3) x SUr(3) symmetry is preserved. In fact, the total 
effective action is invariant under the infinitesimal global chiral transformation 

S F (x,y) -> (l + i^0")S F (x,y)(l-i¥-9 a ), (22) 

S F (x,y) - (l + i^e a )S F (x,y)(l + i^9 a ) (23) 

except for the quark bare mass term and the flavor mixing term which breaks the Ua(1) 
symmetry. 



4 Schwinger— Dyson Equation 

The Schwinger-Dyson equation is derived by the stability condition of the CJT effective action 

6 I\M = o. (24) 
SS Fmn (x,y) 

The detailed procedure is same as in Ref . |T6[| . Introducing the regularized propagators S F (q) := 
f(q 2 )S F (q) and S$(q) := — m ), the SD equation in momentum space becomes 

C r 1 
J{—q ) Jp J{—p ) 

^ P ,fc f(-V 2 )f(-k 2 ) 



x 



W /(-p 2 )/(-^ 2 ) 

tr (DC) [5|? 9l/l ( P )]tr( DC )[^ 2/2 (fc)] (25) 



where the indices m,n, ••• are combined indices m := (a,i,f), n := (b, j, g), ■ ■ ■ with Dirac 
indices a,b, - ■ ■ and color indices z, j, • • • and flavor indices f,g,-- ■■ This equation is shown 
diagramatically in Fig.|l| Generally the quark propagator is parametrised by 

s " fe) = iTO (26) 
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where the index h denotes the flavor. After the Wick rotation, we obtain Ah(—q%) = 1. Then 
the resulting SD equation reads 

g d nMs) r*» irw(r) tb,£) a dlw(l) m-i) (27) 



8vr 4 Jo KJ r + B 2 (-r)Jo I + B 2 (-l) 
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3C F f A lv , _ 2 , N B s (-r) f„, , , r. 



B a (-s) = m s + -— drg (s,r) — r- ^0(r - s) + - r) 

G D Nlw{s) ( A uv , , , rfl-(-r) [ A uv J7 .,. Z5,(-Z) 

^rw(r) ^- ; / dlw(l) - ■ !\. y ; . , (2* 



8tt 4 Jo w r + 52(-r)7o W Z + E 9 2 (-/) 

with s = q E . These integral equations are solvable numerically. 

Since the improved ladder approximation (ILA) model reproduces the asymptotic behavior 
of QCD, the quark mass function can be renormalized so that the solution of the SD equation is 
matched with the QCD quark mass function in the aymptotic region. We renormalize the quark 



mass function properly in the manner described in Ref.[|T6[. The renormalization constants Z m 
and Z ma defined by m q = Z^m q n and m s = Z^m S R are determind by the condition 

M =1 M d M =1. (29) 

om qR m qR =o om sR m gR =0 

Note that the flavor mixing interaction (FMI) does not disturb the asymptotic behavior of the 
ILA model and QCD because of the Gaussian type weight function. 

5 Bethe— Salpeter Equation 

The homogeneous Bethe-Salpeter (BS) equation is derived by 



S 2 T[S_ 



,,v, f -r-i\Xn'm'(y',x';P B ) = (30) 

where 

X n>m>(y',x';P B ) = (0\T^ n ,(y')^ m ,(x')\P) (31) 

denotes the BS amplitude. The normalization condition is (Pb\P'b) = (27t) 3 2Pbo^ 3 (-Pb — 



P'b) an d Pb '■= (y M B + P b ,Pb) is the on-shell momentum of the meson. Introducing the 
regularized BS amplitude by 

XnmWi Pb) ■= f(-q 2 +)XnM P B )f(-d), (32) 

9 
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Figure 2: The diagram for the BS equation. 



(33) 



the BS equation in momentum space becomes 



~iC F I r , ,,w, ln2 , f(ql,kl)iD^(q-k)(^ X R (k-,PBhu)nrn 



k f(-kl)f(-kl) 

+ 2iG n e 9hf 'e fh ' 9 '5 - I - 

+ D Jt Lf(-P 2 )fH 

xtr( DC )[5|?,, /l (p)]{(75) ; , a tr( DC )[ 75X ^,(A ; ; P B )\ + l ba tr^[ X f, f ,(k; P B ))} 



(34) 



which is shown diagramatically in Fig.|2|. For the pseudoscalar state \Pb), the last term in the 
braces does not contribute. 

For the pion, the BS amplitude can be written in terms of four scalar amplitudes as in 



Ref. 16 



X * m (k; P) = 6. 
1 



(a«; 



gf 



j i 



^s(k;P)+4> P (k;PW 

- ba 



P)f 



(35) 



where A" denotes the flavor structure of the pion state. For example, the neutral pion is given 
by a = 3 

/l 0\ 

A 3 = -1 . (36) 
\0 0/ 

On the other hand, for the T] and rj' mesons, the BS amplitudes are written in terms of eight 
scalar amplitudes, 



xit m (k;P) = 5 ji 



(k;P) + ^(k;P)^ + ^(k;P)r 
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75 



ba 



+ 5 Jt 



)gf 
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+ -^(k-P){fk:-mH 



ba 



(37) 



(38) 



where the flavor matrices A", X ss are defined by 
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(39) 



This is because the qq and ss components of the BS amplitude are mixed through the last 
term of the right hand side in Eq . ([34]) which is FMI. The BS equations for the r\ and rj' mesons 
are same. The ground state solution is identified with the r] meson and the first excited state 
solution is identified with the rj' meson. Substituting Eq.(^) or Eq.(^) into Eq. (p4|) , we obtain 
the coupled integral equations. The explicit form is rather complicated and we do not show it. 
Formally the equations can be written as 



Jk 



(40) 



(pA or (ps denotes the set of amplitudes, (fis, (ftp, <Pq, 4>t for the pion and 01?, 0p, 0q , 4>t > 
4> s i , 4>p , 4>n , 4>t f° r t ne V an d ff meson. Instead of solving Eq.( [40|) directly, we solve an eigen- 
value problem 



XMq;Pb)= I M M3 (q,k;P B )Mk;PB) 

Jk 



(41) 



for a fixed P B = M B = —Pbe = ~M E < 0. Then we plot and extrapolate the eigenvalue A as 
a function of P B and search for the on-shell point A = 1. 



6 Decay Constant 



To obtain the decay constant, we need the normalization of the BS amplitude which is derived 



from the inhomogeneous BS equation. In Ref.[|16|] the normalization condition in the momentum 
space is given by 



lim i 



Xnimife PB)Xm 2 n 2 ( ( l', Pb) 







-2P 2 dPv 
11 



cR-l 



S^- 1 (a 

LJ Fmim,2 Vy 



P, 

2' 



1 (42) 



where the integral region is 

I={q\-(q±P B ) 2 < A 2 VY }. 

Using the normalized BS amplitude, the decay constant is obtained by 
1 r 1 



lim 



So f(-C 



-tr 



P^P, P 2 J, S[-qi)S(-qi) 

+(-!(-<£) + f (-<&)*) 
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n n,mm 
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+ | /C n'n,W (_ k + p^ q __._ q __ 
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2GD(75)oi% / w _L_ tr <OC, [S H i/2( , )] 

. (, _ Lf - k 2 - P) - w( -q 2 - -P 2 -k 2 -l 2 + Pl)\ 
2 4 J 
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X£ // 2 /3 e «3 tr ( DC )[(^ i ) 93/3 (0] 
P 



+ | u ,(_( 5+ £)2_ jb 2_^_ u; ^_ 5 2_3 p2 _ A;2 _ z2 _ jPZ 



Xe M e «3 tr ( DC )[(^ f \ /3 (!)] 



(43) 



(44) 



(45) 



(46) 



The term given in terms of E a (q;P) and F a (q;P) represents the correction to the Noether 
axialvector current due to the momentum dependencies of the effective lagrangian. For the 
pion, we can choose the flavor structure matrix \ a of the Noether current in the above formula 
to match that of the BS amplitude. Then we obtain the decay constant f n in the usual sense. 
On the other hand, for the i] or rj' mesons, the flavor structure of the BS amplitude depends on 
the relative and total momenta in general. Therefore we can not fix A" from the flavor structure 
of the BS amplitude. Instead we only consider the decay constants associated with the octet 
(a = 8) and singlet (a = 0) axialvector currents for the r] or r/mesons, i.e., /§ , /q , /g and /q . 
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The fact that the flavor structure of the r]-T]' meson BS amplitudes depend on the relative and 
total momenta means that one cannot define the rj-rj' mixing angle. It can be defined only in 
the limit of neglecting these momentum dependences. 



7 Nambu— Goldstone Solution 

A remark is given here about the Nambu-Goldstone solution. In the chiral limit, the effective 
action is invariant under the 577^(3) x SUr(3) x C/y(l) transformation. Under the dynamical 
breakdown of this symmetry to S77y(3) x Uv(l), we expect eight Nambu-Goldstone (NG) 



solutions. This can be proved by the same procedure as in Ref . [53] . However, we show here the 
existence of these NG solutions directly from the SD equation (^) and the BS equation ([34]). 
Multiplying the 7 5 A a /2 from left and / j 5 X a /2 from right to Eq. (|25|) , we obtain 

A a A a 
{^75 y , S* _1 (g)}nm - {275 — , S* _1 (g)} n m 

x l k f { -/)f { ^) w{ - q2 - p2 - k2 ^ {DC) ^ (P)^ {DC) l S ^fM- (47) 

In the chiral limit the second term of the left hand side vanishes. Furthermore if X a is an octet 
matrix i.e. tr[A°] = 0, then a relation 

[ {^gg' q'hh r /j?l.<72 4. e 9hhJ'9l92 i^)fl \ X s 

I 2 -ret 2 / 1 

= -4 e 9hh e f9i 9 2 6gih ( ^%Il (48) 

leads us to 
A a 

{il^,S§' 1 {q)} nm 

= ~7F?) / P 7F^)^ 2(g| ' p|) ^ (p " g) (v^^yb.) 

+ 2 i G D 5 Jt 5 ba e 9hh e f ^ f 1 , 2 M ~g 2 " P 2 - k ") 

J P ,k f(-p 2 )f(-k 2 ) 

x tr (DC) [^ ffl/l (p)]tr( DC )[ T5 {i T5 ^, Sfrh)}^]. (49) 
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Comparing this with Eq. (p4|) , we obtain the Nambu-Goldstone solution for a = 1, ... ,8 



xL(q; P = 0)= iV{z 7s y , f(-q 2 )S$(q)} (50) 

where N is a normalization constant. It can be also shown that N equals to 1/J'b where fs 
is a decay constant defined in the previous section in the framework of the present ladder like 
approximation. It should also be noted that Eq.(^) does not hold for a singlet A a=0 . 



8 Anomalous PCAC Relation 

The matrix element of the PCAC relation (|17D between a meson state (P\ and the vacuum |0) 
in the ladder like approximation becomes 

-f«M 2 B = 2{m £«]+6 a0 A B (51) 

where f% is the decay constant f n , f% , /g > fs or fo ■ [ m o^i?] an d Ab are defined by 

A. - pltaj, I /( _ 4) 1 /( _e) trfeHfa; PBM(9; p)1 ' (53) 



4™(«;P) := -3G I ,( 75 )^, J ^ | /( _ F) 1 /( _ p) tr' DC »[g« j/2 (t)] 
< - (9 - y) 2 - * 2 " f 2 )^ / '^ s 'tr< TC )[(S?4) 1 ,, / ,(()] 



2 / j L v , 2 



(54) 



respectively. This relation is obtained systematically using the method of Sec. 3 in Ref.j23j. Of 
course, we can also obtain this relation fl5"I|) directly from the SD equation and the BS 
equation (34l). If we employ the BS amplitude in the chiral limit like Eq . ([50"D , it holds that 

[m Q E%\ = m q (qq) Q /f B = m qR (qq) R /f B . (55) 

Eq.([5i~]) with Eq.(^) leads us to the Gell-Mann, Oakes and Renner mass formula 

M 2 B f B ~ -2m q (qq) for a ^ 0. (56) 
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For later use, we define the ratio 

7? (p 2 ) = fB( P E) P E 

which is to be unity at the on-mass-shell point of the Bethe-Salpeter solution. This condition 
is useful in checking the numerical extrapolation procedure. 



9 Numerical Results 

In the present model, there are seven input parameters. Five of them are the parameters of 
the improved ladder approximation (ILA) model of QCD: the current quark mass m q R for 
the up and down quarks, the scale parameter of QCD Aq CD , the infrared cut-off tpp for the 
running coupling constant, the smoothness parameter to and the ultraviolet cut-off Auv- We 
take the value of to from the result of Ref.[jnj, namely, to = —3. As explained there this 
smoothness parameter i is introduced just for the stability of the numerical calculations and 
has no physical meanings. We take Auv = 100 [GeV] because the physical observables depend 
on it rather weakly after the renormalization as far as we use a reasonably large value of Auv- 
The renormalization point /i is taken as \x = 2 [GeV]. The infrared cut-off tip controls the 
strength of the running coupling constant in the low q 2 E region. Therefore its value is directly 
related to the size of the dynamical chiral symmetry breaking. We take tip = —0.5. In the case 
of no FMI, tjp = -0.5 gives -(^V)ij 3 = 259 [MeV] with A QC d = 600 [MeV] in the chiral limit. 

We choose A QCD = 600 [MeV]. Although this value is somewhat larger than the "standard" 
value Aqcd = 100 ~ 300 [MeV], it is necessary for the rainbow gluon self energy to generate 
the dynamical chiral symmetry breaking strongly O, [T(J. Of course if the dynamical chiral 



symmetry breaking is caused by the flavor mixing interaction mainly, we may choose a smaller 
Aqcd- But in such a situation, a careful analysis is necessary. Since the phase transition is of 
the first order, there may exist multiple SD solutions. We will report such results elsewhere 



24J. Therefore in this paper we concentrate only on the case that the chiral symmetry breaking 
is generated mainly by the gluon exchange interaction. 

We have two new parameters associated with the flavor mixing interaction (FMI), i.e., Gd 
and k. Instead of Gd, we use the parameter Iq defined by 

Gd [GeV -5 ] =: — (Ig [GeV -1 ]) 5 . (58) 
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This parameter is chosen freely so that we study the effects of the {7a (1) anomaly on the 77-77' 
system. The k parameter is taken as 

k = 0.7 [GeV~ 2 ]. (59) 

This value corresponds to the form factor of the instanton of the average size p, about 1/3 [fm]. 
The instanton form factor, 
1 r 2 

oc 1 - ±f + . . . (60) 



Xe + P P 

can be identified with the Fourier transformation of the weight function 

F.T. W (g|) = Cexp^-^ cx 1 - || + - - - (61) 

with 

Ak = p 2 . (62) 

The values of the model parameters we use throughout this article are Auv — 100 [GeV], 
Aqcd = 600 [MeV], t = -3, *if = -0.5, p = 2 [GeV] and k = 0.7 [GeV~ 2 ]. 

Let us now discuss the solutions of the SD equation. Our numerical results are shown in 
Table UJ and Figs|J and |4j. As can be seen from them, the chiral symmetry breaking is induced 
mainly by the gluon exchange interaction, and the effect of FMI to the chiral quark condensate 
seems very small. When Iq increases from zero to 2.4, B q (0) increases only 4-6%, and (qq)n 
changes by about 10 to 20 %. One may wonder whether the range of variation is too small for 
its effects to be seen, but, as we will see later, this Iq gives a large mass to 77 and rj'. Since the 
perturbative quark mass contribution to the quark condensate is subtracted in our definition, 
the absolute value of (ss)r is smaller than that of (qq)u (q = u,d). Our results in m q n = 5 [MeV] 
and m sR = 100 [MeV] are (ss)r/ (qq)n = 0.85 and 0.63 for I G = and 2.4 [GeV -1 ] respectively. 
These are in reasonable agreement with the QCD sum rule results: (ss)r/ (uu)r = 0.8 ±0.1 J2"5| 
and (ss)r/ (uu)r = 0.6 ± 0.1 pEfl . The absolute value of the u,d-quark condensate increases as 
Iq increases, while that of the s-quark condensate decreases as I G increases. It is an interesting 
feature in the present model. In the case of the Nambu-Jona-Lasinio (NJL) model with FMI, 
both condensates increase. We do not find out the intuitive explanation of the difference of 
this behavior between the NJL model and the present model. 
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Figure 3: q\ dependences of the solutions of the SD equation in the SU(3) chiral limit with 
I G = 0.0,1.6,2.0,2.4 [GeV- 1 ]. 




Figure 4: q E dependences of the solutions of the SD equation in the SU(2) chiral limit with 
m sR = 100 [MeV] and I G = 0.0, 1.6,2.0,2.4 [GeV" 1 ]. 
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m qR [MeV] 


m sR [MeV] 


ZctGeV- 1 ] 


B q (0) 


B s (0) 


(ftMGeV 3 ] 


{ss) R [GeV 3 ] 








0.0 


1.10 


1.10 


-(0.259) 3 


-(0.259) 3 








1.6 


1.11 


1.11 


-(0.260) 3 


-(0.260) 3 








2.0 


1.11 


1.11 


-(0.262) 3 


-(0.262) 3 








2.4 


1.14 


1.14 


-(0.268) 3 


-(0.268) 3 





100 


0.0 


1.10 


1.22 


-(0.259) 3 


-(0.245) 3 





100 


1.6 


1.11 


1.22 


-(0.261) 3 


-(0.244) 3 





100 


2.0 


1.12 


1.22 


-(0.266) 3 


-(0.241) 3 





100 


2.4 


1.16 


1.24 


-(0.276) 3 


-(0.235) 3 


5 


100 


0.0 


1.11 


1.22 


-(0.259) 3 


-(0.245) 3 


5 


100 


1.6 


1.11 


1.22 


-(0.261) 3 


-(0.244) 3 


5 


100 


2.0 


1.13 


1.23 


-(0.266) 3 


-(0.241) 3 


5 


100 


2.4 


1.17 


1.24 


-(0.276) 3 


-(0.236) 3 



Table 1: Dependences of the solutions of the SD equations at q\ — and the quark condensates 
on the strength of FMI with three sets of the quark masses, (i) m qR = m sR = 0, (ii) m qR = 0, 
m sR = 100 [MeV] and (hi) m qR = 5 [MeV], m sR = 100 [MeV]. Here q represents the u and d 
quarks. 
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2.4 
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101 




5 


100 


0.0 


159 


88 


1.06 


5 


100 


1.6 


158 


90 


1.06 


5 


100 


2.0 


157 


94 


1.05 


5 


100 


2.4 


152 


103 


1.05 



Table 2: Dependences of the solutions of the pion BS equation on the strength of FMI with 
three sets of the quark masses, (i) m q R = m sR = 0, (ii) m qR = 0, m sR = 100 [MeV] and (iii) 
m qR = 5 [MeV], m sR = 100 [MeV]. 

Let us now turn to the discussion of the solutions of the BS equation. Our numerical results 
for the pion are summarized in Table |2|. Here we have not performed the precise parameter 
fittings so as to reproduce the observed pion mass and decay constant since solving the BS 
equation of the non-local interaction requires the rather large computer resouces. We observe 
that the pion mass and decay constant are not sensitive to the flavor mixing interaction. The 
deviation of the slope M%/m qR from the slope —2(qq) R /f% derived from Gell-Mann-Oakes- 
Renner (GMOR) relation is about 16% at m qR = 5 [MeV] and I G = 2.4 [GeV" 1 ]. We consider 
this amount of the deviation of the GMOR relation may come from the Euclid — > Minkowski 
extrapolation, since the ratio TZ defined in Eq.(|5"7|) deviates from unity by 5%, which indicates 
the size of the numerical error for M% and /„- in the extrapolation procedure. 

The BS solutions for the rj and rj' mesons are given in Tables^ and |j. Since the BS equation 
is homogeneous, the absolute sign of the BS amplitudes, and therefore the decay constants, 
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cannot be determined. We choose the sign of f$ (/o) to be positive for 77 (77'). The masses of 77 
and 7/ and their decay constants depend strongly on the flavor mixing interaction. Especially, 
the 7/ meson mass seems sensitive to the flavor mixing interaction. This is in contrast to the 
pion result. {7a(1) symmetry breaking gives a large effect on the 77 and 77' sector. 

In order to see the effects of the flavor mixing, we introduce the mixing angles for the 77 and 
77' mesons, 

—pf = tan9 v , -j-r = tan 6^/. (63) 
J& Jo 

The results are presented in Tables ^| and [|. Since the flavaor structure of the 77-7/ meson 
BS amplitudes depend on the relative and total momenta, the above definitions of the mixing 
angles are the kinds of the averaged quantities. 

In the SU(3) symmetry limit, no flavor mixing occurs and 9 Tj = 9 V > = 0. On the other 
hand, in the broken SU(3) case without FMI, 77 and 77' are in the ideally mixed states, i.e., 
9 = arctan(- v / 2) = -54.7°. The mixing angles 9 V and 9 v i increase as FMI becomes strong. In 
the case of m q R = 5 [MeV], m s # = 100 [MeV] and Iq = 2.4 [GeV -1 ], the reasonable values of 
the M v and M v i are obtained, the calculated M v is 7% smaller than the observed M v and the 
calculated M v i is 11% larger than the observed M v i. In this model parameters, the calculated 
mixing angle for the 77' meson 9 V > is more than 5/3 times of the calculated mixing angle for the 
77 meson 9 r] . It means that the momentum dependences of the flavor structures of the 77 and 77' 
mesons are not so small and the momentum independent treatment of the 77-77' mixing angle is 
rather questionable. 

The last column of Tables 0, |3] and |4] gives the ratio 1Z in Eq.(^) for finite quark mass. As 
it should be 1 at the on-mass-shell point identically, it is a good indicator of the ambiguity, 
or error, coming from the extrapolation from the Euclidean kinematics to the Minkowski on- 
mass-shell kinematics. Here we carry out the quadratic extrapolations of the eigenvalue A 
and the ratio TZ and the linear extrapolation of the decay constant / from the Euclid region 
to the on-mass-shell point. In the case of the weak FMI, our extapolation procedure works 
rather well. The quality of the extrapolation can be seen in FigJ^. For heavier meson masses, 
1Z deviates from 1 significantly. This indicates an extrapolation error. In fact, for Mb > 700 
[MeV] the extrapolation becomes very difficult in the quadratic extrapolation. For instance, TZq 
for Iq = 2.4 [GeV -1 ] shown in Table |3] largely deviates from 1. Fig.|6] shows the extrapolation in 
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this case, where the lines of A and 1Z% are almost straight but the curve of TZq is not. This may 
be a reason why TZq deviates from unity. We have performed the extrapolation of TZq by using 
a rational function which is shown in Fig|| by the dotted line and the result is improved well. 
As for the rj' meson, the extrapolation of the eigenvalue, the rations and TZq are shown in 
Fig.0. 

From Table [| one can see that in the chiral limit, rf state is a pure flavor singlet state and 
has finite mass. This mass is due to FMI and it means 7/ is not the Goldstone boson. One of the 
interesting question is that how much rf loses the the Goldstone boson nature. In the present 
range of the Ua{^) breaking interaction strength, the flavor singlet pseudoscalar meson state 
has mass from 194 [MeV] to 634 [MeV]. On the other hand the decay constant changes only 
less than 8%. Further studies of the 7/ properties such as the decay amplitudes are necessary 
in order to understand the nature of the rf meson. 

We plot the rf meson mass as a function of the L^(l) breaking parameter Iq in Fig.§. The 
effect of the mixing of u, d quark component seems to be negligible and the 77' mass grows 
rapidly from I G ~ 2.0 [GeV^ 1 ]. 

The 77-meson properties have been studied extensively in the three-flavor NJL model with 
the instanton induced Ua(1) breaking flavor mixing interaction (FMI) in Ref.|ll| and it has 
been shown that the r/-meson mass, the rj —>■ 77, 77 — > 7, [i~fi + and rj —>■ 7r°77 decay widths 
are reproduced well with the rather strong FMI. In this case the contribution from FMI to 
the dynamical mass of the up and down quarks is about 44% of that from the usual C/l(3) X 
Ur(3) invariant four-quark interaction. In contrast with it, the contribution from FMI to the 
dynamical quark mass is very small in the present study. To make the situation clear, let us 
compare the strength of FMI in the present case with that in the NJL model case. The naive 
way is to compare the following two quantities, 

(64) 







in the present model and 

pA 167r 4 







dq ^ G °> (65) 



in the NJL model. In this manner Iq = 2.4 [GeV corresponds to G e § = 0.73, which is rather 
close to the value determined in Ref.|ll[, Gff = 0.7. It means the strength of FMI is almost 
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fV 

Jo 
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[MeV] 


[MeV] 


[GeV- 1 ] 


[MeV] 


[MeV] 


[MeV] 


[deg] 












0.0 





86 








— 











1.6 





87 


















2.0 





89 


















2.4 





95 















100 


0.0 





50 


70 


-54.7 









100 


1.6 


203 


56 


68 


-50.5 


1.03 


1.04 





100 


2.0 


351 


71 


58 


-39.2 


1.00 


1.03 





100 


2.4 


495 


106 


22 


-11.7 


1.00 


0.62 (1.07) 


5 


100 


0.0 


159 


51 


72 


-54.7 


1.05 


1.05 


5 


100 


1.6 


258 


57 


69 


-50.4 


1.02 


1.03 


5 


100 


2.0 


387 


75 


58 


-37.7 


1.00 


1.05 


5 


100 


2.4 


511 


113 


19 


-9.5 


1.01 


0.55 (1.00) 



Table 3: Dependences of the ??-meson solutions of the coupled channel BS equation on the 
strength of FMI with three sets of the quark masses, (i) m qR = m sR = 0, (ii) m qR = 0, 
m sR = 100 [MeV] and (iii) m qR = 5 [MeV], m sR = 100 [MeV]. The values in the brace are 
obtained by the rational extrapolation. 
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m qR 


m sR 


Ig 


M v 


fV 
J8 


Jo 


9 V 


Til 




[MeV] 


[MeV] 


[GeV- 1 ] 


[MeV] 


[MeV] 


[MeV] 


[deg] 












0.0 








86 

















1.6 


194 





87 







1.04 








2.0 


350 





88 







1.02 








2.4 


634 





94 







1.06 





100 


0.0 


723 


-101 


72 


-54.7 


1.05 


1.05 





100 


1.6 


731 


-99 


76 


-52.5 


1.05 


1.07 





100 


2.0 


764 


-88 


88 


-45.0 


1.06 


1.10 





100 


2.4 


1020 


-41 


118 


-19.2 


1.09 


1.30(1.32) 


5 


100 


0.0 


723 


-101 


72 


-54.7 


1.05 


1.05 


5 


100 


1.6 


732 


-98 


77 


-51.8 


1.05 


1.07 


5 


100 


2.0 


777 


-86 


90 


-43.7 


1.06 


1.11 


5 


100 


2.4 


1060 


-36 


122 


-16.4 


1.09 


1.33 



Table 4: Dependences of the 7/-meson solutions of the coupled channel BS equation on the 
strength of FMI with three sets of the quark masses, (i) m qR = m sR = 0, (ii) m qR = 0, 
m sR = 100 [MeV] and (iii) m qR = 5 [MeV], m sR = 100 [MeV]. The value in the brace is 
obtained by the rational extrapolation. 
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Eigen Value and Ratios 




Euclid Mass Square [GeV 2 ] 

Figure 5: Extrapolation of the eigenvalue and the ratios for the r\ meson with m q n = 0,m s ^ = 
100 [MeV] and I G = 1.6 [GeV- 1 ]. 



Eigen Value and Ratios 




Euclid Mass Square [GeV ] 



Figure 6: Extrapolation of the eigenvalue and the ratios for the t] meson with 
m sR = 100 [MeV] and I G = 2.4 [GeV 1 ]. 
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Figure 7: Extrapolation of the eigenvalue and the ratios for the meson with m q R = 0, 
m sR = 100 [MeV] and I G = 2.4 [GeV" 1 ]. 
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Figure 8: Iq dependence of the mass of the rj' meson with m q n = and m s ^ = 100 [MeV]. 



same in the both cases. It is not clear why the contribution from FMI to the dynamical quark 
mass is so different in the two cases. 

It should be noted here that as can be seen from Fig.^], FMI operates on the pseudoscalar 
meson states as the effective 4-quark interaction reduced by contracting a quark-antiquark pair 
into a quark condensate. Therefore the effective strength of FMI is not Gd but Gd (<?<?)• In our 
present model parameters, the dynamical chiral symmetry breaking is mostly driven by the one- 
gluon exchange type interaction. If one reduces the strength of the one-gluon type interaction 
in the infrared region, the quark condensate becomes smaller and the effective strength of FMI 
on the meson states becomes weaker. Therefore there is a possibility of taking rather strong 
FMI without destroying the success of the present description of the i] and rj meson masses. 
In the case where FMI is dominant in the infrared region, we expect that the chiral phase 
transition becomes the first order and there may exist multiple solutions of the SD equation. 
In such a situation, more careful analyses are required, which will be reported elsewhere |24 . 
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10 Summary and Conclusions 



The improved ladder approximation of QCD has successfully described the low-energy proper- 
ties of QCD |13], plj. In this article, we have studied the rj and rj' mesons in this approach. 
It is expected that the Ua(1) anomaly plays an important role in the rj and rj' mesons and 
therefore we have introduced the instanton induced Ua(1) breaking 6-quark determinant in- 
teraction [^, [7]] in the improved ladder approximation model of QCD. We have derived the 
Schwinger-Dyson (SD) equations for the light quark propagators and the Bethe-Salpeter (BS) 
equations for the pion, r\ and rf in the lowest order (rainbow-ladder) approximation using the 
Cornwall- Jackiw-Tomboulis (CJT) effective action formulation [plf| . 



Using the same model parameters of the running coupling constant used in [I6|, we have 
obtained reasonable values of M n , M v , M v /, f n and (qq)R with a relatively weak flavor mixing 
interaction (FMI), for which the chiral symmetry breaking is dominantly induced by the soft- 
gluon exchange interaction. It is in contrast with the Nambu-Jona-Lasinio (NJL) model results, 
where about 1/3 of the dynaminacl quark mass is due to FMI |iT |. 

As far as we know, the rf BS equation which includes the running coupling aspect of QCD 
and the effect of the Z7a(1) anomaly has not been solved so far. In the case of the NJL model, 
the rf mass has unphysical large imaginary part associated with the unphysical decay channel 
rf — > qq. On the other hand, the present model predicts a bound rj' although it may not perfectly 
confine quarks. The bound state is obtained as the quark mass function becomes rather large, 
B(—q 2 E ) ~ 1 [GeV] at g# = 0, in this model, while it is independent of the momentum in the 
NJL model. 

Since the flavor structure of the rj-rj' meson BS amplitudes depend on the relative and total 
momenta, one cannot define the rj-rj' mixing angle unambiguously. It can be defined only in the 
limit of neglecting these momentum dependences and it should be examined whether such an 
approximation is reasonable. Our numerical results indicate that the momentum independent 
treatment of the flavor mixing is rather questionable. 

In the chiral limit we can define the rjo decay constant /o without any ambiguity. (Here 770 
means the pure flavor singlet state.) Our numerical results show that (|/o| — f n )/ fn — 0.012 in 
our parameter range though M m changes from 194 [MeV] to 634 [MeV]. There is no low-energy 
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theorem for the r] decay constant /o since the Ua(X) symmetry is explicitly broken by the 
anomaly. Therefore the present result of f should contain the information of the low-energy 
dynamics of QCD. 

The present result is our first step towards quantitative understanding of the flavor mixing 
interaction. With the BS solutions in hand, we may calculate static properties and decay 
amplitudes of rj and rj'. Our model is regarded as a low energy effective theory, which is 
consistent with chiral symmetry, its spontaneous breakdown and the Ua(X) anomaly. It should 
be stressed that the approximation used in solving and renormalizing the amplitudes also 
respect these symmetry properties. Thus our approach is suitable for further studies of the 77 
and i]' systems, which are desirable in order to clarify the role of the ?7a(1) anomaly in the 
low-energy QCD. 
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